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A Dalitz plot analysis of approximately 12500 D° events reconstructed in the hadronic decay 
D° —> K K + K~ is presented. This analysis is based on a data sample of 91.5 fb _1 collected 
with the BABAR detector at the PEP-II asymmetric-energy e + e~ storage rings at SLAC running at 
center-of-mass energies on and 40 MeV below the T(4S) resonance. The events are selected from 
e + e~ — > cc annihilations using the decay D* + — > D°ty + . The following ratio of branching fractions 
has been obtained: 

BR = -i _ ; = (15.8 ± 0.1 (stat.) ± 0.5 (syst.)) x 10" 2 . 

Estimates of fractions and phases for resonant and non-resonant contributions to the Dalitz plot are 
also presented. The ao(980) — » KK projection has been extracted with little background. A search 
for CP asymmetries on the Dalitz plot has been performed. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 
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I. INTRODUCTION 

The Dalitz plot analysis is the most complete method 
of studying the dynamics of three-body charm decays. 
These decays are expected toproceed through interme- 
diate quasi-two-body modes |l| and experimentally this 
is the observed pattern. Dalitz plot analyses can pro- 
vide new information on the resonances that contribute 
to observed three-body final states. 

In addition, since the intermediate quasi-two-body 
modes are dominated by light quark meson resonances, 
new information on light meson spectroscopy can be ob- 
tained. Also, old puzzles related to the parameters and 
the internal structure of several light mesons can receive 
new experimental input. 

Puzzles still remain in light meson spectroscopy. There 
are new claims for the existence of broad states close to 
threshold such as n (800) and cr(500) @. The new evi- 
dence has reopened discussion of the composition of the 
ground state J PC = ++ nonet, and of the possibility 
that states such as the ao(980) or / (980) may be 4- 
quark states due to their proximity to the KK thresh- 
old 3] . This hypothesis can only be tested through an ac- 
curate measurement of branching fractions and couplings 
to different final states. In addition, comparison between 
the production of these states in decays of differently fla- 
vored charmed mesons D°(cu), D + {cd) and D+(cs) can 
yield new information on their possible quark composi- 
tion. Another benefit of studying charm decays is that, 
in some cases, partial wave analyses are able to isolate 
the scalar contribution almost background free. 

This paper focuses on the study of the three-body D° 
meson decay 

where the K° is detected via the decay Kg — > 7r + 7r~ . All 
references in this paper to an explicit decay mode, unless 
otherwise specified, imply the use of the charge conjugate 
decay also. 

This paper is organized as follows. Section II briefly 
describes the BABAR detector, while Section III gives de- 
tails on the event reconstruction. Section IV is devoted 
to the evaluation of the efficiency and the measurement 
of the branching fraction is reported in Section V. Sec- 
tion VII deals with a partial wave analysis of the K + K~ 
system, while sections VI, VIII, IX and X describe the 
Dalitz plot analysis. 



tAlso with Universita della Basilicata, Potenza, Italy 
J Deceased 



II. THE BABAR DETECTOR AND DATASET 

The data sample used in this analysis consists of 91.5 
ftT 1 recorded with the BABAR detector at the SLAC 
PEP-II storage rings. The PEP-II facility operates nom- 
inally at the T(AS) resonance, providing collisions of 
9.0 GeV electrons on 3.1 GeV positrons. The data 
set includes 82 fb _1 collected in this configuration (on- 
resonance) and 9.6 fb~ collected at a cm. energy 40 
MeV below the T(45) resonance (off-resonance). 

The following is a brief summary of the components 
important to this analysis. A more complete overview of 
the BABAR detector can be found elsewhere Q. The in- 
teraction point is surrounded by a five-layer double-sided 
silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) filled with a gas mixture of helium and isobutane, 
all within a 1.5-T superconducting solcnoidal magnet. In 
addition to providing precise spatial hits for tracking, 
the SVT and DCH measure specific energy loss dE/dx, 
which provides particle identification for low-momentum 
charged particles. At higher momenta (p > 0.7 GeV/c) 
pions and kaons are identified by Cherenkov radiation 
observed in the DIRC, a detector designed to measure 
Cherenkov angles of photons internally reflected in the 
radiator. The typical separation between pions and kaons 
varies from 8a at 2 GeV/c to 2.5a at 4 GeV/c. 

III. EVENT SELECTION AND D° 
RECONSTRUCTION 

This analysis includes a measurement of the branching 
ratio: 

ER= T(D°^K°K+K-) 
T{D° -> K°Tr+ir-) ' 

Therefore the selection of both the data samples corre- 
sponding to 

D°->K\ + 7[-, (1) 

and 

D° -» K°K+K- (2) 

is described. The two final states are referred to collec- 
tively as K°h+h-. 

The decay D* + — ► D°ir + is used to distinguish between 
D° and D° and to reduce background. For example, the 
Cabibbo-favoured decays under study are 

D*+ -> D°n+ 

-> K°ir+iT- 



D*- -> D°ir- 

-> K°TT + TT- 

The charge of the slow ir^ from D* decay (referred to 
as the slow pion ng) identifies the flavor of the D° and 
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K° (for the latter, ignoring the small contribution from 
doubly-Cabibbo-suppressed decay of the D°). 

A D° — > K s h + h~ candidate is reconstructed from 
a K s — > 7T + 7r _ candidate plus two additional charged 
tracks, each with at least 12 hits in the DCH. The slow 
pion is required to have momentum less than 0.6 GeV/c 
and to have at least 6 hits in the SVT. In addition, all 
the tracks are required to have transverse momentum 
Pt > 100 MeV/c and, except for the Kg decay pions, 
to point back to the nominal collision axis within 1.5 cm 
transverse to this axis and within ± 3 cm of the nominal 
interaction point along this axis. 

A K s candidate is reconstructed by means of a ver- 
tex fit to a pair of oppositely-charged tracks with the K° 
mass constraint. The reconstructed Kg candidate is then 
fit to a common vertex with all remaining combinations 
of pairs of oppositely-charged tracks to form a D° candi- 
date vertex. K s candidates are further required to have 
a flight distance greater than 0.4 cm with respect to the 
candidate D° vertex. The D° candidate is then combined 
with each slow pion candidate, and fit to a common D* 
vertex, which is constrained to be located in the interac- 
tion region. In all cases the fit probability is required to 
be greater than 0.1%. 

To reduce combinatorial background, a D° candidate 
is required to have a center of mass momentum greater 
than 2.2 GeV/c. 

Kaon identification is performed by combining dE/dx 
information from the tracking detectors with associated 
Cherenkov angle and photon information from the DIRC. 
The resulting efficiency is above 95% for kaons with less 
than 3 GeV/c momentum that reach the DIRC. 

Each D° sample is characterized by the distributions 
of two variables, the invariant mass of the candidate D°, 
and the difference in invariant mass of the D* + and D° 
candidates 

Am = m(K°h + h-Ti+) - m{K°h + h-). 

The distributions of Am for those candidates for which 
m(K°h + h~) is within two standard deviations of the 
D° mass value are shown in Fig. Qfa) and Fig. ^c) 
for K°tt + it~ and K a K + K" respectively. Strong D* 
signals are apparent. Fits to these distributions pro- 
duce consistent means and widths for the two channels: 
Am = 145.41 ±0.01 MeV/c 2 , a = 304±4KeV/c 2 in the 
case of decay channel K°ir + ir~ (statistical errors only). 

The K°h + h~ mass distributions for candidates that 
fall within ± 600 KeV/c 2 (two standard deviations) of 
the central value of the Am distribution are shown in 
Fig. □» and Fig. QJd) for K°tt+tt- and K°K+K~ re- 
spectively. 

Fitting the K°h + h~ mass distributions using a linear 
background and a Gaussian function for the signal gives 
the following mass and width (statistical errors only) for 
the decay D° -> Wn+Tr-: 

m = 1863.65 ± 0.06 MeV/c 2 ; a = 6.10 ± 0.02 MeV/c 2 

and for D° ->K°K + K--. 




Am(K' r 7T + tv~) (GeV/c 2 ) 
(c) 




m(K° n + tT) (GeV/c 2 ) 
I (d) 



0.1452 0.149 

Am(K° K* K~) (GeV/c 2 ) 




1.81 1.86 1.91 

m(K° K* K") (GeV/c 2 ) 



K°h+h 



FIG. 1: a) and c) The Am distributions for D° 
candidates, for events in which the K°h + h~ invariant mass 
is within two standard deviations of the D° mass value. The 
arrows indicate the region of Am used to select the D° can- 
didates, b) and d) K°h + h~ mass distributions for events in 
which Am is within 600 KeV/c 2 of the mean Am value for 
signal events. The arrows indicate the region of m(K°h + h~) 
used to produce the Am distributions. 



m = 1864.74 ±0.03 MeV/c 2 ; a = 3.37 ±0.03 MeV/c 2 . 

The mass resolution for reaction (2) is much better than 
that for reaction (1) because of the much smaller Q- value 
involved (380 MeV/c 2 compared to 1088 MeV/c 2 ). The 
w 1 MeV/c 2 shift between the two mass measurements 
is within the expected systematic error and is due to the 
different kinematics of the two D° decay modes. 



IV. EFFICIENCY 

The selection efficiency for each of the D° decay modes 
is determined from a sample of Monte Carlo events in 
which each decay mode was generated according to phase 
space (i.e. such that the Dalitz plot is uniformly popu- 
lated). These events were passed through a full detector 
simulation based on the GEANT4 toolkit [f| and sub- 
jected to the same reconstruction and event selection pro- 
cedure as were the data. The distribution of the selected 
events in each Dalitz plot is then used to determine the 
relevant reconstruction efficiency. Typical Monte Carlo 
samples used to compute these efficiencies consist of 2 
x 10 5 generated events. Each Dalitz plot is divided into 
small cells and the efficiency distribution fit to a third- 
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order polynomial in two dimensions. Cells with fewer 
than 100 generated events were ignored in the fit. The 
resulting \ 2 per degree of freedom (x 2 /NDF) is typically 
1.1 using ps 500 cells. The fitted efficiencies are shown 
m Fig. Using the weighting procedure described in 
the next section the weighted efficiencies values (17.94 ± 
0.25)% for K°Tr+n- and (16.56 ± 0.38)% for K°K+K- 
are obtained. The above errors include the uncertainties 
on the weighting procedure. 

V. BRANCHING FRACTIONS 

Since the two K°h + h~ decay channels have similar 
topologies, the ratio of branching fractions, calculated 
relative to the K°ir + ir~ decay mode, is expected to have 
a reduced systematic uncertainty. This ratio is evaluated 
as 

2-hb,V e (x,y) 

where Ni(x, y) represents the number of events measured 
for channel i, and Ci(x,y) is the corresponding efficiency 
in a given Dalitz plot cell (x,y) . 

To obtain the yields and measure the relative branch- 
ing fractions, each K°h + h~ mass distribution is fit 
assuming a double Gaussian signal and linear back- 
ground where all the parameters are floated, as shown 
in Fig. njb) and Fig. IHd). The number of signal events 
is calculated as the difference between the total num- 
ber of events from the fit and the integrated linear back- 
ground function in the same mass range. The region used 
is within ±6<r of the D° mass. Selecting events within 
three standard deviations of the central value of the Am 



distribution, the fits give the following yields for the two 
channels. 

D° -> K°n+TT- : N = 92935 ± 305 

D° -^K°K+K- : N = 13536 ± 116. 

Systematic errors take into account effects due to the 
use of selection regions for the Am distribution, the use of 
particle identification, the different fitting models used to 
subtract the background, Kg reconstruction, and uncer- 
tainties in the calculation of the efficiency on the Dalitz 
plot due to Monte Carlo statistics. The resulting sys- 
tematic error is dominated by the uncertainty due to ef- 
ficiency correction on the Dalitz plot. 

The resulting ratio is: 

br _ T(D°^K°K+K-) _ 

(15.8 ± 0.1 (stat.) ± 0.5 (syst.)) x 10~ 2 

to be compared with the PDG value of: (17.2 ± 1.4) x 
10 -2 |6J. The best previous measurement of this branch- 
ing fraction comes from the CLEO experiment (136 
events for reaction (2)), which obtains the value BR = 
(17.0 ±2.2) x 10~ 2 jf. 

The branching ratio measurements have been validated 
using a fully inclusive e + e _ — * cc Monte Carlo simulation 
incorporating all known D° decay modes. The Monte 
Carlo events were subjected to the same reconstruction, 
event selection and analysis procedures as for the data. 
The results were found to be consistent, within statistical 
uncertainty, with the branching fraction values used in 
the Monte Carlo generation. 
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FIG. 3: Dalitz plot of D° -> K°K + K~ 



VI. DALITZ PLOT FOR D 



o 



K°K + K~ 



Selecting events within ±2<r of the fitted D° mass 
value, a signal fraction of 97.3% is obtained for the 
12540 events selected. The Dalitz plot for these D° — > 
if°iir + Ar _ candidates is shown in Fig. 03 In the if+if - 
threshold region, a strong 0(1020) signal is observed, to- 
gether with a rather broad structure. A large asymmetry 
with respect to the K°K + axis can also be seen in the 
vicinity of the 0(1020) signal, which is most probably the 
result of interference between S and P-wave amplitude 
contributions to the K + K~ system. The /o(980) and 
ao(980) S- wave resonances are, in fact, just below the 
K + K~ threshold, and might be expected to contribute 
in the vicinity of 0(1020). An accumulation of events 
due to a charged ao(980) + can be observed on the lower 
right edge of the Dalitz plot. This contribution, how- 
ever, does not overlap with the 0(1020) region and this 
allows the K + K~ scalar and vector components to be 
separated using a partial wave analysis in the low mass 
K + K~ region. 



VII. PARTIAL WAVE ANALYSIS 

It is assumed that near threshold the production of the 
K + K~ system can be described in terms of the diagram 
shown in Fig.^J The helicity angle, 9k, is then defined as 



0=0 





FIG. 4: The kinematics describing the production of the 
K + K~ system in the threshold region. 



the angle between the K + for D° (or K~ for D°) in the 
K + K~ rest frame and the K + K~ direction in the D° (or 
K°) rest frame. The K + K~ mass distribution has been 
modified by weighting each D° candidate by the spherical 
harmonic Y® (cos 6k) (L=0-4) divided by its (Dalitz-plot- 
dependent) fitted efficiency. The resulting distributions 
(Y L ) are shown in Fig. and are proportional to the 
K + K~ mass-dependent harmonic moments. It is found 
that all the \Y®) moments are small or consistent with 
zero, except for (r o °), (Yf) and (F 2 °). 

In order to interpret these distributions a simple partial 
wave analysis has been performed, involving only S- and 
P-wave amplitudes. This results in the following set of 
equations 



47r(y °) = S 2 + P 2 



2|S||P|cos0 S p 



(3) 



v^(r 2 °) = -|p 2 , 

where 5* and P are proportional to the size of the S- and 
P-wave contributions and 4>sp is their relative phase. 
Under these assumptions, the (Y§) moment is propor- 
tional to P 2 so that it is natural that the 0(1020) ap- 
pears free of background, as is observed. This distribu- 
tion has been fit using the following relativistic P-wave 
Breit-Wigncr. 

For a resonance r — > AB, BW(m) is written as 



BW(m) 



Fr 



mf. 



m 



AB 



iT AB m r 



(4) 
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FIG. 5: The unnormalized spherical harmonic moments (Yj? ) as functions of K + K invariant mass. The histograms represent 
the result of the full Dalitz plot analysis. 



where for a spin J=l particle F r is the Blatt-Weisskopf 
damping factor [j| 



F r = 



y/i + (Rq AB ) 2 

and R has been fixed to R = 1.5 GeV -1 . In Eq. (4): 

q r m AB 

where qab {Qr) is the momentum of either daughter in 
the AB (r) rest frame. 

The fit yields the following parameters: 

m = 1019.63 ± 0.07, I> = 4.28 ± 0.13 MeV/c 2 

in agreement with PDG values (statistical errors only). 
The fit is shown in Fig. |SJ 

A strong S — P interference is evidenced by the rapid 
motion of the (Yx) moment in Fig. 0m the </>(1020) mass 
region. 



The above system of equations (3) can be solved di- 
rectly for S 2 , P 2 and cos 4>SP- However, since these am- 
plitudes are defined in a D° decay, it is necessary to cor- 
rect for phase space. This has been achieved by using 
the K + K~ and K°K + mass spectra obtained from the 
Monte Carlo generation of D° decays to K°K + K~ ac- 
cording to phase space. The D° mass distribution has 
been generated in this Monte Carlo as a Gaussian hav- 
ing the experimental values of mass and mass resolution. 

The phase space corrected spectra are shown in Fig.0 

The distributions have been fitted using the following 
model: 

• The P-wave is entirely due to the cf>(1020) meson 
(Fig. Eta)). 

• The scalar contribution in the K + K~ mass projec- 
tion is entirely due to the a (980)° (Fig.I2»). 

• The K°K + mass distribution is entirely due to 
a (980)+ (Fig. He)). 
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FIG. 6: (Y®) spherical harmonic moment as a function of the 
K + K~ effective mass. The line is the result from the fit with 
a relativistic spin-1 Breit Wigner. 



• The angle 4>sp (Fig.Efd)) is obtained fitting the S, 
P waves and cosfisp with c ao BW ao + c^BW^e 101 . 



Here BW ao and BW<p are the Breit- Wigner describ- 
ing the ao(980) and 0(1020) resonances. 

The do (980) scalar resonance has a mass very close to 
the KK threshold and decays mostly to Tjir. It has been 
described by a coupled channel Breit Wigner of the form: 



BW ch (a )(m) 



9kk 



il-m 2 - iip^g 2 ^ + Pkk9 2 Rk ) 



(5) 



where p(m) = 2q/m while g^ and g R x describe the 
ao(980) couplings to the r/ir and KK systems respec- 
tively. 

The best measurements of the an (98 0) p arameters 
come from the Crystal Barrel experiment 10], in pp an- 
nihilations, and are the following: 



m =999±2 MeV/c 2 , 5?)7r =324±15 (MeV) 
%-=1.03±0.14. 

9kk 



1/2 



This corresponds to a value oig RK = 329_±27 (MeV) 1 / 2 . 

Since in the current analysis only the KK projections 
are available, it is not possible to measure mo and g v „. 
Therefore, these two quantities have been fixed to the 
Crystal Barrel measurements. The parameter g RR , on 
the other hand, has been left free in the fit. The result 
is (statistical error only): 

9RK = 464 ± 29 (MeV) 1 / 2 . 



Figure [TJe) shows the residual ao(980) phase, obtained 
by first computing <fisp in the range (0,7r) and then sub- 
tracting the known phase motion due to the 0(1020) 
resonance. The fit gives a value of a relative phase 
a = 2.12 ± 0.04 and has a \ 2 /NDF=167 /92. The fit is 
of rather poor quality, indicating an undetermined source 
of systematic uncertainty comparable with the statistical 
uncertainty. However the issue related to the determina- 
tion of g RK will be rediscussed in the complete Dalitz 
plot analysis described in section VIII. 

The entire procedure has been tested with Monte Carlo 
simulations with different input values of the ao(980) pa- 
rameters. The partial wave analysis performed on these 
simulated data yielded the input value of g RKl within 
the errors. 

In this fit the possible presence of an /o(980) contri- 
bution has not been considered. This assumption can 
be tested by comparing the K + K~ and K°K + phase 
space corrected mass distributions. Since the /o(980) 
has isospin 0, it cannot decay to K°K + . Therefore an ex- 
cess in the K + K ~ mass spectrum with respect to K°K + 
would indicate the presence of an /o(980) contribution. 

Figure HO compares the K + K~ and K°K + mass distri- 
butions, normalised to the same area between 0.992 and 
1.05 GeV/c 2 and corrected for phase space. It is possible 
to observe that the two distributions show a good agree- 
ment, supporting the argument that the Jo (980) contri- 
bution is small. Notice that the enhanced K°K + signal 
level above 1.1 GeV/c 2 is the result of the 0(1020) reflec- 
tion. 

The resulting scalar components of the K + K~ and 
K°K + mass distributions, corrected for phase space, are 
tabulated as a function of mass in Table [I] 



VIII. DALITZ PLOT ANALYSIS OF 

D° -> K°K + K- . 

An unbinncd maximum likelihood fit has been per- 
formed for the decay D° — > K°K + K~ in order to use 
the distribution of events in the Dalitz plot to determine 
the relative amplitudes and phases of intermediate reso- 
nant and non-resonant states. 

The likelihood function has been written in the follow- 
ing way: 



(3'G(m)e(ml,rriy) 



Ey C i C *j A i A *j 



I T,i,j CiC*AiA*e{m 2 x , m 2 y )dmldm 2 y 



+ (1-/3) (6) 

In this expression, (3 represents the fraction of signal ob- 
tained from the fit to the mass spectrum and e(m 2 , m 2 ) 
is the fitted efficiency on the Dalitz plot. The Gaussian 
function G(m) describes the D° line shape normalised 
within the ±2cr cutoff used to perform the Dalitz plot 
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FIG. 7: Results from the K + K~ Partial Wave Analysis corrected for phase space, (a) P-wave strength, (b) S-wave strength, 
(c) m(K° K + ) distribution, (d) cos^sp in the 0(1020) region, (e) cj>sp in the threshold region after having subtracted the fitted 
0(1020) phase motion shown in (d). The lines correspond to the fit described in the text. 



analysis. It is assumed that the background events, de- 
scribed by the second term in Eq. (6), uniformly populate 
the Dalitz plot. This assumption has been verified by ex- 
amining events in the D° side bands. The output from 
the fit is the set of complex coefficients c, . 

In Eq. (6), the integrals have been computed using 
Monte Carlo events while taking into account the effi- 
ciency on the Dalitz plot. The branching fraction for the 
resonant or non-resonant contribution i is defined by the 
following expression: 

= \ Cl \ 2 f\Ai\ 2 dm 2 x dm 2 y 

Ej,k c 3 c U A i A l dm l dm l' 

The fractions fa do not necessarily add up to 1 because 
of interference effects among the amplitudes. The errors 
on the fractions have been evaluated by propagating the 
full covariance matrix obtained from the fit. 

The phase of each amplitude is measured with respect 
to A'°a (980)° which gives the largest contribution. The 
amplitudes A, are represented by the product of com- 



plex Breit-Wigner BW(m) (Eq. (4)) and angular terms 
T{Q) [13: 

A = BW{m) x T{Q). 

The /o(980) resonance has been described using a cou- 
pled channel Breit-Wigner function with parameters 
taken from the WA76 [3, E791 0, and BES 13]. The 
ao(980) has been parametrized using the results from the 
partial wave analysis discussed above. 

The parameters of the </>(1020) meson have been fixed 
to the values obtained from the fit to the (F 2 °) moment 
described earlier. The non-resonant contribution (NR) is 
represented by a constant term with a free phase. 

Systematic errors on the fitted fractions have been 
evaluated by making different assumptions in the fits. 
For example, in one test, the efficiency on the Dalitz plot 
has been set to a constant value. In other tests the reso- 
nance parameters of /o(980), ao(980) and /o(1400) have 
been fixed to values obtained from a variety of experi- 
ments. 
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TABLE I: K + K~ and K°K+ scalar mass projections cor- 
rected for phase space in arbitrary units. 



800 




0.95 1 1.05 1.1 1.15 

m(K K) (GeV/c 2 ) 

FIG. 8: Comparison between the phase-space-corrected 
K + K~ and K°K + normalised to the same area in the mass 
region between 0.992 and 1.05 GeV/c 2 . 



The doubly-Cabibbo-suppressed contribution (DCS) 
A" + a (980)~, whose presence should appear like an 
ao(980) _ in the wrong sign combination K°K~ , has been 
also included in the ht. 



IX. RESULTS FROM THE DALITZ PLOT 
ANALYSIS. 

The D° — > K°K + K~ Dalitz plot projections together 
with the fit results are shown in Fig. [§] 

Figure [7| shows the fit projections onto the (Vj) mo- 
ments. The fit produces a reasonable representation of 
the data for all of the projections. The \ 2 computed on 
the Dalitz plot gives a value of \ 2 /NDF=983/77<i. The 
sum of the fractions is 130.7 ± 2.2 ± 8.4%. The regions of 
higher \ 2 are distributed rather uniformly on the Dalitz 
plot. Attempts to improve the fit quality by including 
other contributions did not give better results. One par- 
ticular problem found in these fits is that including too 
many scalar amplitudes caused the fit to diverge, produc- 
ing a sum of fractions well above 200% along with small 
improvements of the fit quality. 

The final fit results showing fractions, amplitudes and 
phases are summarised in Table ITU For K°j b(980) and 
K + a (980)~ (DCS), being consistent with zero, only the 
fractions have been tabulated. 

The results from the Dalitz plot analysis can be sum- 
marised as follows: 



mass / c J 


K+K 




: \ i\ 


n noo 
0.988 


644 ± 105 




n not 
0.992 


474 


± 


52 


0(0 ± to4 


n aaa 
0.99b 


417 


± 


37 


AO A _L OO 

484 ± 82 


i nnn 
t.000 


392 


± 


37 


4t4 ± DO 


1.004 


304 


± 


35 


noo _i_ a o 
282 ± 48 


t.OOo 


299 


± 


33 


OO 1 _L ACi 

33t ± 49 


1.012 


259 


± 


39 


010 _l_ o o 

2t3 ± 38 


t.Oto 


240 


± 


62 


OOE A- OO 

2oo ± 38 


1.020 


178 


± 


84 


ion _i_ o o 

189 ± 33 


t.024 


210 


± 


45 


ICO A- OO 

to3 ± lo 


t.028 


178 


± 


30 


1 (17 A- OO 

t9/ ±61 


t.032 


157 


± 


23 


1 on A- OK 

iz9 ± 2o 


1.036 


164 


± 


19 


140 ± 25 


1.040 


147 


± 


20 


102 ± 21 


1.044 


135 


± 


17 


117 ± 22 


1.048 


139 


± 


15 


132 ± 23 


1.052 


126 


± 


13 


114 ± 22 


1.056 


101 


± 


14 


119 ± 22 


1.060 


101 


± 


12 


108 ± 20 


1.064 


104 


± 


11 


72 ± 17 


1.068 


120 


± 


12 


63 ± 15 



• The decay is dominated by D° K°a {980) , 
D° -> K°cj)(1020) and D° -> ^"^(980)+. 

• The /o(980) contribution is consistent with zero, 
even after assuming various / (980) lineshape pa- 
rameters [13 E3 

• The DCS contribution is consistent with zero, re- 
gardless of the an (980) parametrisation. 

• The remaining contribution is not consistent with 
being uniform, but can be described by the tail 
of a broad resonance, for example the Jo (1400) 
which peaks well outside the phase space. It is 
not possible to derive its parameters from our data, 
but several parametrizations have been tried, in 
particular those from J/ip decays 13] and from 
Df — > 7r + 7r + 7r~ getting in all cases improved 
fits. 

• In one of the fits the /o(1400) contribution has been 
replaced by a non-resonant contribution, obtaining 
a fraction of 2.5 ± 0.5%. However the likelihood 
value for this fit was worse A(21ogL) = 56. 

For the if°/o(980) and DCS contributions upper limits 
have been computed. Combining statistical and system- 
atic errors in quadrature, the following 95 % c.l. upper 
limits on the fractions have been obtained: 

BF(D° -> i?°/o(980)H K+K-)) < 2.1% 



BF(D° -» K+a (980)-(-^ K°R-)){DCS) < 2.5%. 
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TABLE II: Results from the Dalitz plot analysis of D° — » K°K + K . The fits have been performed using the value of K = 464 
(MeV) 1 / 2 resulting from the partial wave analysis. 



Final state 


Amplitude 


Phase (radians) 


Fraction (%) 


F"a (980) u 


1. 


0. 


66.4 ± 1.6 ± 7.0 


A"V(1020) 


0.437 ± 0.006 ± 0.060 


1.91 ± 0.02 ± 0.10 


45.9 ± 0.7 ± 0.7 


A"-a (980)+ 


0.460 ± 0.017 ± 0.056 


3.59 ± 0.05 ± 0.20 


13.4 ± 1.1 ± 3.7 


Jf u /o(1400) 


0.435 ± 0.033 ± 0.162 


-2.63 ± 0.10 ± .71 


3.8 ± 0.7 ± 2.3 


if u /o(980) 






0.4 ± 0.2 ± 0.8 


K+a (980)- 






0.8 ± 0.3 ± 0.8 


Sum 






130.7 ± 2.2 ± 8.4 
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A test has been performed by leaving g^ K as a free 
parameter in the Dalitz plot analysis. In this test the 
other parameters describing the a (980) (m and g m ) 
have been allowed to vary within their measurement er- 
rors from the Crystal Barrel experiment. The resulting 
central value of g^K 1S 473 (MeV) 1 / 2 with a maximum 
deviation of 39 (MeV) 1 / 2 , in good agreement with the 
value obtained using the partial wave analysis. The dif- 
ference between the values, added in quadrature with the 
above maximum deviation, has been taken as an estimate 
of the systematic error: 



9RK = 473 ± 29 (stat.) ± 40 (syst.)(MeV) 



1/2 



m 2 (K° K + ) (GeV 2 /c 4 



This value differs significantly from the Crystal Barrel 
measurement. An improvement of this measurement can 
be foreseen by adding data from the a (980) — * i]ir decay 
mode such as D° — > K^rjir . This D° decay mode has 
been studied by the CLEO Q experiment (with rather 
limited statistics) finding a D° — > i\r o a o (980)° dominant 
contribution. 

A large uncertainty is included in the upper limit on 
the presence of Jo (980) in this D° decay mode due to the 
poor knowledge of the /o(980) parameters. A small signal 
of /o(980) is indeed present (in this case as a shoulder) 
in the D° -> K°tt + tt- as shown in Fig. [TU1 

Dalitz plot analyses of this D° decay channel have been 
performed by BABAR [H andJBelle [l| finding (w 5.5%) 
as decay fraction for D° — > A"°/ o (980). However, a reli- 
able estimate of the expected contribution of the /o (980) 
in D° — > K°K + K~ decay is not possible until more ac- 
curate measurements of the /o(980) parameters and cou- 
plings become available. This can be performed, for ex- 
ample, by using high statistics samples of Df — » KKtt + 
and Z?+ — * 7r + 7r + 7r~ decays. 



X. SEARCH FOR CP ASYMMETRIES ON THE 
DALITZ PLOT. 



FIG. 9: Dalitz plot projections for D° -> K°K+K~. The 
data are represented with error bars; the histogram is the 
projection of the fit described in the text. 



A search for CP asymmetries on the Dalitz plot has 
been performed. Table IIIII shows the results from the 
Dalitz plot analysis performed separately for D° and D a . 
Notice that in these two fits good values of x 2 /NDF have 
been obtained. 
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TABLE III: Results from the Dalitz plot Analysis of D° -> K°K + K~ separated for D° and D°. 



Decay mode 


fraction (%) 


amplitude 


phase (radians) x 2 /NDF 


D" 


-> K°a (980) 


66.5 ± 2.0 


1. 


0. 


671/649 


D° 


-> K°a (980) 


66.3 ± 2.0 


1. 


0. 


643/646 


D u 


-> K>(1020) 


46.3 ± 0.8 


0.438 ± 0.009 


1.93 ± 0.03 




D° 


K>(1020) 


45.6 ± 0.8 


0.435 ± 0.009 


1.88 ± 0.03 




D u 


-» JC-a (980) + 


13.2 ± 1.3 


0.456 ± 0.025 


3.58 ± 0.07 




D° 


— ► K+a (98Q)- 


13.6 ± 1.3 


0.463 ± 0.025 


3.59 ± 0.07 




D u 


^K»f {U00) 


4.1 ± 0.9 


0.451 ± 0.047 


-2.58 ± 0.13 




D° 


^K°fo{UQQ) 


3.6 ± 0.9 


0.421 ± 0.038 


-2.68 ± 0.14 





2500 




m(n + tT) (GeV/c 2 ) 



FIG. 10: 7r + 7r~ effective mass from D° -> K°n + n~. The 
arrow indicates the position of the /o(980). 



We do not observe any statistically significant asym- 
metries in fractions, amplitudes, or phases between _D° 
and D°. 



XI. SUMMARY 

A Dalitz plot analysis of the D° hadronic decay D a — > 
K°K + K~ has been performed. The following ratio of 
branching fractions has been obtained: 



(15.8 ± 0.1 (stat.) ± 0.5 (syst.)) x 10~ 2 . 
The Dalitz plot analysis indicates that the channel is 
dominated by D° -> i?°a (980) , D° -» 7?V(1020) and 
D° -» ^^00(980)+. The a (980) -» KK lineshapc has 
been extracted with little background. 

The Dalitz plot analysis of D° and D° do not show any 
statistically significant asymmetries in fractions, ampli- 
tudes, or phases. 



BR = 



T(D° -> K K+K~) 
T(D° K°ir+TT-) 
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